Pt-incorporated anatase TiO2(001) surface for solar cell applications : First-principles 

density functional theory calculations 
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First-principles density functional theory calculations were carried out to determine the low energy 
geometries of anatase TiC>2(001) with Pt implants in the sublayers as substitutional and interstitial 
impurities as well as on the surface in the form of adsorbates. We investigated the effect of such a 
systematic Pt incorporation in the electronic structure of this surface for isolated and interacting 
impurities with an emphasis on the reduction in the band gap to visible region. Comprehensive 
calculations, for lxl surface, showed that Pt ions at interstitial cavities result in local segregation, 
forming metallic wires inside, while substitution for bulk Ti and adsorption drives four strongly 
dispersed impurity states from valence-bands up in the gap with a narrowing of ~1.5 eV. Hence, 
such a contiguous Pt incorporation drives anatase into infrared regime. Pt substitution for the 
surface Ti, on the other hand, metallizes the surface. Systematic trends for 2x2 surface revealed 
that Pt can be encapsulated inside to form stable structures as a result of strong Pt-O interactions 
as well as the adsorptional and substitutional cases. Dilute impurities considered for 2x2 surface 
models exhibit flat-like defect states driven from the valence bands narrowing the energy gap suitable 
to obtain visible light responsive titania. 

PACS numbers: 68.43.Bc, 68.43. Fg 



I. INTRODUCTION 

The wide-gap semiconductor titania (TiC^) has raised 
great interest primarily because of its catalytically active 
surfaces, long-standing stability, non-toxicity, and avail- 
ability of single-crystal samples**^ 3 - The most common 
phases of titania are known to be anatase, rutile, and 
brookitc, among which anatase phase proves to be the 
most promising for photoelectrochemistry^ visible-light 
photocatalysis^ii rocking-chair lithium batteries^ and 
optoelectronics.— 

Hcngcrcr et a/.— studied the stability of anatase T1O2 
(101) and (001) facets and found that it is possible to ob- 
tain clean and structurally perfect anatase surfaces. Sin- 
gle crystals of anatase TiC>2 exhibit stronger photocat- 
alytic activity than rutile phase titania . 11 ' 12 The natural 
crystallographic (001) surface of anatase is most often 
considered for catalytic applications among its various 
facets^ 

Titania is active under UV irradiation while it is al- 
most inert to solar spectrum by absorbing only 2%-3 % 
of the sun light. Hence, narrowing the band gap to visible 
range is particularly important for practical photocatal- 
ysis. Such applications utilize excess electrons incorpo- 
rated by various impurities. Resultant defect states fall 
in the band gap of TiC-2 sensitizing it for visible light- 
induced catalytic purposes. Usual way of defect forma- 
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tion is achieved by oxygen vacancy which reduces the 
oxide. However, merely the vacancy driven defect states 
are not enough for high level of activity. 

Many attempts have been made to functionalize the 
titania surfaces for solar cell applications by impurities in 
the form of ion doping or dye sensitization! 7 - : 14 : 15 : 16 : 17 : 18 
Moreover, by acting as charge trap sites, such impurities 
are proposed to help in retarding the fast charge-hole 
recombination rates, which inherently exist in most of 
the semiconductors as TiC-2. 

Doping TiC>2 with metallic as well as nonmctallic el- 
ements has been extensively studied for powdered pho- 
tocatalysts. Anatase phase of TiC>2 doped with N, S, 
C, and B has been reported to exhibit relatively high 
level of visible activity i 18 ' 19 : 20 ' 21 : 22 : 23 Supportingly, Wang 
et ai.fii in their theoretical study, showed that N dop- 
ing narrows the band gap of TiC>2 by bringing impurity 
states in the vicinity of valence-band maximum (VBM). 
Such an enhancement can also be obtained by transition- 
metal ion doping which reduces the gap allowing visible- 
light absorption by providing intraband states near the 
conduction-band or valence-band edges . 14 i 24 i 25 Yet, it de- 
pends on the role of dopants as recombination centers or 
as charge traps. For instance, Co 3+ and Al 3+ impurities 
serve as electron-hole recombination sites, significantly 
decreasing the photoreactivity. Pt ion doping has been 
successfully shown to enhance visible activity function- 
ing as charge generation centers which produce free and 
trapped chargesi 7 -^ 

Metal ion doping in the form of Ti substitution has 
been proposed also for titania based dye-sensitized solar 
cells (DSSC) which gain visible light activation through 
dye molecule surfactants. In this manner, DSSC photo- 
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voltaic efficiencies have been found to be remarkably bet- 
ter for doped TiC>2 by preventing injected dye electron re- 
combination between the electrolyte and the substrate.— 
On the other hand, Pt incorporation is not merely limited 
to powdered photocatalytic systems. Recently, Kitano et 
al^L successfully developed visible-light responsive Pt- 
loaded Ti02 thin-film photocatalysts which achieve sep- 
arate H2 and O2 evolution from water without requiring 
dye sensitization. 

Encapsulation of Pt in titania due to the strong metal 
support interaction (SMSI) under reducing-gas atmo- 
sphere has been reported by Pesty et al.*& Later, Zhang 
et al£L showed that neutral Pt atoms can thermally dif- 
fuse into TiC>2 lattice under oxidizing atmosphere. They 
also argue that these diffused Pt atoms can substitute for 
Ti 4+ when oxidized to Pt 2+ or else they form interstitial 
impurities as well. 

In this paper, we have investigated the effect of Pt in- 
corporation in both the lattice and the electronic struc- 
ture of stoichiometric anatase TiC>2(001) in the form of 
strongly interacting and noninteracting impurities on and 
in the surface. The aim is to shift the activity to visible 
region for solar cell operation. This is done so by band 
gap narrowing which is useful for photovoltaic devices. 
We have considered Pt as adsorbates on the surface, and 
as substitutional and interstitial impurities in the sub- 
surface layers encapsulated by the slab. 



II. METHOD 

We performed total-energy and electronic-structure 
calculations using the Vienna Ab Initio Simulation Pack- 
age (VASP) implementation^ of the gradient-corrected 
[Perdew-Burke-Ernzerhof (PBE)] (Refi^) density func- 
tional theory (DFT). The electron-ion interaction has 
been described by the projector augmented waves (PAW) 
metho d 30 i 31 using plane-wave basis sets. 

The naturally occurring rutile and anatase polymorphs 
of titania are basically formed as a result of different mod- 
ifications of the same TiOg unit. This building block is 
arranged as a distorted octahedron with a Ti cation at 
the center and six oxygen at the vertices. The stack- 
ing pattern of these octahedra results in simple tetrag- 
onal (st) conventional unit cells for both of rutile and 
anatase polymorphs. However, in the case of anatase, 
the primitive unit cell is a body-centered tetragonal ( bet) 
Bravais lattice. Since the conventional cell for anatase 
contains two bet units, the calculations assuming an st 
unit cell might cause misleading deductions such as the 
cnergy-band-gap type. We considered bulk properties of 
anatase phase of TiC>2 as the starting point, both for 
building up the surface slab models and for obtaining 
the bulk-projected electronic structures. In this manner, 
we calculated, for example, the bulk lattice parameters 
to be a = 3.801 A, c = 9.468 A, and u = 0.2095 with 
a (IAi/amd) space-group symmetry using bet unit 
cell. These results agree well with the corresponding ex- 



perimental values^ which were reported as a = 3.785 A, 
c = 9.514 A, and u = 0.208. Besides, our PAW-GGA 
calculations for these structural properties are consistent 
with the other available theoretical result a 33 i 34 ' 35 i 36 and 
exhibit a better agreement with the experimental find- 
ings. 

For the stoichiometric anatase TiO2(001) surface we 
considered an oxygen terminated supercell model involv- 
ing six Ti02 layers with a vacuum region of -13 A. These 
zigzag-like Ti02 layers consist of three atomic layers in 
which bridging oxygen atoms are out of the level Ti plane. 
Each of the oxygen at the back surface of such a bulk 
termination needs 1/3 electrons to be saturated. This 
cannot be accomplished by hydrogenation with integral 
charge. Indeed, it leads hydrogen driven surface states 
to appear just above the VBM as a consequence of the 
excess charge induced by this hydrogenation. Therefore, 
instead of saturating the back surface we chose a vir- 
tually symmetrical slab model which produces the same 
electronic properties coming from the top and the bot- 
tom (001) surfaces. We call it as virtually symmetrical 
because it lacks perfect mirror symmetry along the axis 
perpendicular to the surfaces. In this sense it is not a 
trivial surface to model. On the other hand, our choice 
bears no mistakes since we use plane- wave basis sets with 
periodic boundary conditions. 

In order to elucidate the role of Pt incorporation on 
the electronic behavior of anatase TiO2(001) surface, we 
considered strongly interacting and isolated impurities on 
and inside the corresponding supercells. Pt adsorption on 
(lxl) surface, for instance, corresponds to 1 monolayer 
(ML) coverage so that the shortest Pt-Pt distance is at- 
tained. When one considers Pt with (2x2) construction, 
namely, 0.25 ML, Pt-Pt distance becomes 7.53 A on the 
surface which leads to almost isolated impurities. 

Because of the symmetrical nature of the slab model, 
supercell thickness becomes important particularly when 
Pt penetration depth increases. Implanted Pt atoms 
from the top and back surfaces should not interact in- 
side the slab. In addition, the central part of the model 
slab must possess bulk- like properties. In this manner, 
even though none of the atoms were fixed to the bulk 
positions in the geometry optimizations, the atoms at 
the central part retained their original bulk positions for 
shallow enough Pt penetration. Moreover, the number of 
layers have been chosen so that increasing the slab height 
by one more layer did not alter the calculated results sig- 
nificantly. However, our tests showed that six Ti02 layers 
are not enough for Pt impurities placed deeper than 3.4 
A from the surface oxygen which corresponds to the sec- 
ond Ti02 sublayer. Therefore, we have used eight Ti02 
layer thick slab model for such Pt impurities inside the 
slab. 

Our convergence tests showed that the electronic wave 
functions can be expanded into plane waves up to an 
energy cutoff value of 400 eV and that the surface 
Brillouin-zone integrations can be carried out with a 
fc-point sampling of (8x8x1) and (4x4x1) Monkhorst- 
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Pack meshes^! for (lxl) and (2x2) surface unit cells, 
respectively. In both cases, these choices gave a total- 
energy convergence up to a tolerance value smaller than 
0.1 meV. All geometry optimization calculations were 
carried out using conjugate-gradient algorithm based on 
the reduction in the Hellman-Feynman forces on each 
constituent atom to less than 10 meV/A. 



III. RESULTS AND DISCUSSION 

Anatase phase of bulk titania is found to have an in- 
direct gap of 2.08 eV with the valence-band top being 
at about two-thirds of the way along FX. Direct gap is 
slightly larger, 2.12 eV. The calculated gap values are un- 
derestimated as expected due to inherent shortcomings 
associated with the local density approximation (LDA), 
while the experimental gap is 3.20 eV. 38 ' 39 A more recent 
work reports this value to be 3.5 eV4^ The calculational 
underestimation can be overcome by incorporating self- 
energy corrections going beyond the ground-state DFT. 
Indeed, Thulin and Guerra^ 3 - reported a quasi-particlc 
corrected gap of 3.79 eV which is overestimated and 
closer to the result 3.68 eV obtained by Calatayud el al^ 
using nonlocal B3LYP functional. A comparison of our 
bulk band structure (not shown) with the ones reported 
in these references revealed a perfect match except the 
gap width, hence, indicating a scissors type rigid shift of 
the whole conduction-band. Therefore, despite the un- 
derestimated gap, the band structures presented in this 
work can be considered within this context. 
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FIG. 1: (Color on-line) Atomic arrangements of the anatase 
TiO2(001)-2x2 structure. Side views of the (a) ideal and (b) 
relaxed surface unit cells, only up to eight atomic layers are 
shown here. Ti and O atoms are denoted by light gray and 
red (dark) balls, respectively. All bond lengths are given in 
angstroms and all angles are presented in degrees. 



The (001) plane of anatase polymorph is known to 
be catalytically important as it constitutes commercial 
catalystsi 41 ' 42 Yet, a theoretical investigation on the elec- 
tronic structure of this surface is still needed. Therefore, 
we first consider the oxygen terminated (001) clean sur- 
face. The (lxl) unit cell forms in the shape of a square 
with a side of 3.765 A over this plane [see the top view 



in Fig.[3ja)]. In fact, the (2x2) periodicity has also sim- 
ilar formation [see the top view in Fig. Ufa)]. Because 
of the bulk termination, the surface layer is composed of 
fivefold coordinated Ti and twofold coordinated O atoms 
which are denoted as Ti5c and 02c in Fig. [TJ respectively. 
Bridging 02c atoms make two equidistant bonds with 
Ti5c's along [100] possessing a mirror plane symmetry on 
the ideal surface [Fig.rjja)]. The outward position of 02c 
forms a Ti5c-02c-Ti5c angle of 151.7°. In the ideal geom- 
etry, equatorial bond lengths are 1.94 A while the axial 
bonds measure 2.02 A. This conspicuously symmetrical 
structure gets distorted upon relaxation in consistency 
with the previous theoretical studies, 4 ^. The Ti5c-02c 
bond lengths become significantly inequivalent being 2.16 
and 1.79 A mainly because of the relatively larger dis- 
placement of surface oxygen along [100] compared to the 
Til and Ti2. Moreover, Ti5c-02c-Ti5c angle reduces to 
144.4° since 02c's move outward the surface plane while 
the relaxation of Ti5c's in the reverse direction is substan- 
tially larger than the other atoms near the surface. This 
relaxation pattern applies to both of the surface unit cells 
as all calculated values related to the atomic rearrange- 
ments presented in Table U for (lxl) surface are consis- 
tent with those of (2x2). Also, in agreement with the 
results of previous theoretical studies , 43 ' 44 our geometry 
optimizations gave negligibly small displacements along 
[010]. Lazzeri et al.^ additionally, reported that the pla- 
nar 02-Til-03-Ti2 ring becomes slightly skewed making 
a dihedral angle of 6.0° with the (100) plane. However, 
as can readily be seen in Fig. QJb) , our calculations sug- 
gest that the Til-03 bond shortens to 1.96 A and makes 
an angle of 2.46° with (100) plane whereas 02-Ti2 bond 
has a skew angle of 6.0° and a length of 2.00 A. This 
can be attributed to the high reactivity of the surface 
with an indication that not only 02c but also threefold 
coordinated 02 and 03 give contribution to the surface 
electronic properties. 



TABLE I: Geometric structure of anatase TiO2(001) surface. 
Atomic labels refer to Fig. [T] Calculated values for atomic 
displacements with respect to ideal bulk positions are given 
in A. 



Surface 


Label 


[100] 


[010] 


[001] 


lxl 


Ol (02c) 


0.194 


-0.011 


0.034 




02 


0.223 


-0.010 


-0.038 




03 


-0.117 


-0.011 


-0.031 




Til (T5c) 


0.003 


-0.009 


-0.091 




Ti2 


0.011 


-0.011 


-0.011 


2x2 


Ol (02c) 


0.195 


-0.011 


0.034 




02 


0.224 


-0.010 


-0.037 




03 


-0.117 


-0.011 


-0.030 




Til (Ti5c) 


0.004 


-0.010 


-0.088 




Ti2 


0.015 


-0.013 


-0.012 



An ideal-like structure in which atomic relaxations are 
only observed in the direction perpendicular to the sur- 
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face is found to be energetically higher by 0.09 and 0.37 
eV per unit cell for (lxl) and (2x2), respectively. The 
existence of this symmetry preserving structure was also 
reported by Calatayud and Minot^ for (lxl) as well 
as for larger unit cells. Minimum energy structure can- 
not be reached unless the optimization is started from a 
slightly distorted ideal geometry. Relatively small atomic 
displacements as presented in Table U lead to weak en- 
ergy differences between the minimum energy structure 
and the ideal geometry. Therefore, these weak energy 
differences can be seen as an indication of high thermo- 
dynamic stability for anatase (001). Indeed we calculated 
the surface energies in both of the (lxl) and (2x2) cases 
to be 0.92 J/m 2 and 0.94 J/m 2 for relaxed and unrelaxed 
surfaces, respectively. The surface energy is defined as 



1 

2A 



^TiQ 2 



pbulk 
n£ "Ti0 2 



where £/ri0 2 is the total energy of the slab and nE^Q 2 
refers to the energy of the bulk supercell containing an 
equal number of TiC>2 units as the slab. A corresponds 
to the exposed unit-cell area while the factor of 1/2 ap- 
pears because the slab is symmetrical having two faces. 
Our calculated values are in good agreement with the 
GGA-PBE results (££^=0.98 and E££ el =1.12 J/m 2 ) 
of Lazzeri et alM and, with the GGA-PW91 results 
(££f rf= o.89 and £ s u u n / f el =1.04 J/m 2 ) of Calatayud and 
Minot^L. LDA results are reported to give systemati- 
cally larger surface energiesi 43 i 45 ' 46 The lack of accurate 
experimental measurements to compare with prevents us 
to discuss which exchange-correlation scheme does more 
reliably describe the surface properties. 



(a) 




(b) 
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FIG. 2: (Color on-line) Energy bands for the clean anatase 
TiC>2(001) for (a) lxl and (b) 2x2 surfaces. Some of the im- 
portant LDOS contributions to TDOS are shown. Projection 
of bulk continuum is also depicted as shaded areas. 

The major surface bands derived from the valence- 
bands are spilled out into the energy-band gap. Figure [5] 
shows the energy bands for the clean surface of anatase 
TiO2(001), having (lxl) and (2x2) periodicities, and 
the bulk band continuum (shaded regions) projected on 
the corresponding surface Brillouin zones, respectively. 
Fermi level is at 0.40 eV relative to the valence-band top 



and the surface states are filled causing the gap to nar- 
row down to 1.68 eV (from the bulk gap value of 2.08 
eV). These results for (lxl) surface are in consistency 
with those of (2x2) as presented in Fig. [2{b) and in 
Table [III The gap is indirect with surface state hav- 
ing a maximum at M point for (lxl) while it is di- 
rect with the corresponding maximum at T for (2x2) 
due to the surface Brillouin zone folding. Conduction 
band makes a minimum at T in both cases and is bulk- 
like with all the surface solutions being resonance states 
within the bulk band region except for the pocket-states. 
The right-most panel shows the total density of states 
(TDOS) and important local contributions [local density 
of states (LDOS)] to it. The surface states in the gap de- 
rived mainly from the twofold coordinated surface oxygen 
(02c) are seen to have the form of a jump discontinuity 
at Ep, typical of two-dimensional van Hove singularity 
in the DOS due to the critical point at M (or K). The 
contribution of 02 to the LDOS, being comparably much 
less than that of 02c, comes from the lower part of the 
surface states closer to the valence-band top. The DOS 
from 03 (and similarly from Ti5c) extends to the valence 
bulk continuum affecting the surface states even less than 
02 docs. Evidently, the DOS analysis at the upper region 
of the valence-bands indicates that the energy levels of 
02c lying higher than those of fully coordinated oxygen 
must show higher reactivity. 

We have also calculated the (lxl) and (2x2) clean sur- 
face band structures for the ideal-like geometries in which 
Ti5c-02c bond symmetry is preserved. In this case, the 
surface states move upward in the gap region reducing 
the energy-band gap by 0.29 eV. Moreover, 02c driven 
states elevate from the bulk valence-band up into the gap 
increasing the number of available surface states. These 
elevations are not always rigid. For instance, the highest 
lying surface state exhibits some differences between the 
ideal-like and the minimum-energy configurations. At T 
point, this surface state coincides with the valence-band 
top in the lowest energy structure as shown in Fig.OJa), 
whereas it lies 0.29 eV above the valence-band in the 
case of the ideal-like configuration, attesting a promi- 
nent change in the character of this state. Therefore, 
these differences altogether suggest that the relaxations 
of atoms near the surface, noticeably, influence the sur- 
face electronic band structure, in contrary to what was 
asserted by a previous theoretical work4£ 

In practice, titania surface is covered with less than a 
monolayer or even dilute metal adatom concentration in 
catalysis applications. Hence, we systematically studied 
Pt implantation in anatase TiO2(001) with (lxl) and 
extensionally with (2x2) periodicities starting from the 
surface layer up to seven atomic sublayers corresponding 
to a depth of ~5.1 A. First, for Pt doped Ti02 surfaces, 
we substitute Pt atoms for the fivefold Ti atoms (Til), 
closest to the surface layer, and in place of the bulk Ti 
atom (Ti2) at the second Ti02 layer. We refer to these 
substitutional cases as (si) and (s2), respectively. Sec- 
ond, we considered Pt atoms positioned at interstitial 
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sites in oxygen atomic layers where Pt is strongly coordi- 
nated with the nearest-neighbor oxygen. The trend from 
strongly interacting Pt impurities at 1 ML to isolated 
ones at 0.25 ML coverages is expected to help predict 
the more dilute experimental situations. 



A. Contiguous impurities 

Pt-implanted TiO2(001)-l x 1 system has been consid- 
ered with all possible adsorptional, substitutional, and 
interstitial configurations starting from above the surface 
to inside the slab. Pt impurities in and on the lxl surface 
are separated from each other by 3.76 A. Although this 
is much larger than the Pt dimer length, it still main- 
tains a distance close enough for a strong interaction. 
Pt is found to be stable at two adsorption cases which 
are presented in Fig. [3l In the first case, Pt adsorbate 
binds to both Ti5c and 02c at bridge position forming 
a Ti5c-Pt-02c angle of 51°. Pt-Ti5c bond length is 2.55 
A which is significantly close to Pt dimer length while 
Pt-02c distance is 1.97 A. Structural properties of Pt 
for all the cases studied in this work such as Pt-0 and 
Pt-Ti bond lengths as well as the Pt penetration depths 
are presented in Table [TTJ The twofold binding of Pt 
to Ti5c and 02c compensates their undercoordination 
by charge transfer. Hence, the system undergoes a re- 
laxation in favor of reducing the difference between the 
two inequivalent bond lengths (Ti5c-02c) on the surface 
layer. (See Fig.EJa)). Similarly both of the Til-03 and 
Ti2-02 bonds get equal in length to a value of 2.01 A 
making a dihedral angle, with the (100) plane, of 5.9° 
and of 4.6°, respectively. Ti5c-02c-Ti5c angle, on the 
other hand, increases slightly to 146.1° due mainly to 
the downward relaxation of 02c relative to its nearest- 
neighbor Ti5c atoms. In fact, the surface does not show 
a significant reconstruction upon Pt adsorption, for this 
case, and is characterized dominantly by the relaxation 
of surface oxygen. 

In the second case as shown in [Fig. [3Jb)] Pt is four- 
fold coordinated with the surface atoms. It makes two 
equidistant bonds with 02c's along [010] being 1.96 A 
in length while Pt-Ti bond distances are 2.71 and 2.61 
A. This small difference in bond lengths stems from 
the stronger coordination of Pt with neighboring 02c's. 
This strong interaction forms O-Pt-0 lines along [010] 
as shown in the top view of Fig. [3jb). In this geometry, 
Pt is 0.22 A closer to the surface relative to the previ- 
ous case. Ti5c-02c bonds are slightly inequivalent being 
2.21 and 2.10 A similar to the case of the clean surface. 
Nevertheless, Ti5c-02c-Ti5c angle gets subtly reduced to 
121.4° prevailingly as a result of the elevation of 02c in 
a strong coordination with the Pt adsorbate. The struc- 
ture in case (b) is energetically 0.47 eV more favorable 
than that of case (a) because of the increased coordina- 
tion of Pt with the surface atoms. Pt binding energy 
(BE) is calculated to be 2.70 and 2.93 cV/atom for the 
cases (a) and (b), respectively. 




FIG. 3: (Color on-line) Pt on anatase TiO 2 (001)-l x 1 surface. 
Front, side, and top views for the two adsorption cases: (a) 
on the bridging oxygen bond and (b) off the bridging oxygen 
bond. The surface unit cell is indicated in (a) right. 



When Pt is substituted for Ti5c, which is referred as si, 
Pt-03 and Ti2-02 axial bond lengths increase substan- 
tially to 2.18 and 2.09 A, respectively. Besides, Pt-02c 
equatorial bonds become almost equal being 1.95 and 
1.93 A with Pt-02c-Pt angle which reads 152.5°. In this 
structure, Pt-03-Ti2-02 side ring as a whole is skewed 
making a 5.2° dihedral angle with (100) plane. Pt is per- 
fectly aligned with nearest-neighbor oxygen in the same 
lanes over the surface along [100] and [001] directions 
while being at the different atomic layers. 

On the other hand, Pt substitution for bulk Ti (Ti2), 
namely, s2, results in a structure which reflects simi- 
lar topological characteristics with the ideal clean slab 
[shown in Fig. HJa)]. Pt-Ti2 replacement at lxl unit 
cell drives the structure from the relaxed to an ideal-like 
geometry. Even though Pt replaces Ti2, the axial bond 
lengths extend slightly to 2.02 A. Similarly, the difference 
in bond distances of Pt with the neighboring fully coordi- 
nated oxygen become negligibly small being 1.95 A. The 
only substantial displacement with respect to the relaxed 
structure is obtained for 02c which moves up a little bit 
so that Ti5c-02c bond lengths become symmetrized with 
a value of 1.96 A forming an isosceles Ti5c-02c-Ti5c tri- 
angle having an obtuse angle of 146.9°. 

Pt can also be considered at the interstitial sites in 
between the fully coordinated level oxygen. Starting 
from such a configuration, the surface expansively recon- 
structs with a local segregation as a result not only of 
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FIG. 4: (Color on-line) In the upper row: (si) energy bands 
for Pt impurity atom substituted for surface-Ti ion in the 
anatase TiO2(001)-lx 1 surface and (s2) those for Pt impu- 
rity atom substituted for the second-layer Ti ion. Some of 
the important LDOS contributions to TDOS are shown. Pro- 
jection of bulk continuum is also depicted as shaded areas. 
In the lower row: (a) energy bands due to interstitial Pt on 
anatase TiO2(001)-l x 1 surface as in Fig. E^a) and, (b) that 
as in Fig. 13b). 



the SMSI, but also of the stress induced by the web of 
closely spaced Pt impurities inside the slab. For instance, 
when we place Pt in between two 02 atoms (referring to 
Fig. [J), it pushes the first Ti0 2 layer (Ol-Til-02 group) 
upward breaking Til-03 and Ti2-02 bonds. Meanwhile, 
Pt moves in between 01 and 03 forming new equidis- 
tant axial bonds (1.98 A) on a straight line which makes 
an angle of 26.6° with the (100) plane. Since 01 and 02 
coordination numbers interchange, 02 is further elevated 
up forming zigzag Ti5c-02c pattern, this time along [010] 
direction with bond lengths of 1.85 and 2.08 A. In sum- 
mary, Pt interstitial at 02 level for (lxl) surface seg- 
regates the first and the second Ti02 layers by 2.0 A 
relative to the separation of those in the clean surface. 

The coordination of Pt with nearest-neighbor oxygen 
in these model cases signifies the strength of Pt-0 elec- 
trostatic interaction so that Pt is very dominant in dis- 
turbing the surface atomic arrangements. These impu- 
rity driven rearrangements modify the electronic band 
structure to a significant extent by bringing new defect 
states as well as perturbing the already existing surface 
states that originate from the distorted lattice bondings. 

The electronic structure for the geometries in Figs. 
3(a) and 3(b) with Pt ions as surface impurities are pre- 
sented in Figs. 4(a) and 4(b), respectively. Moreover, in 



Fig. 4(sl) and 4(s2), those for the substitutional impu- 
rity at the surface (replacing Til) and that in the bulk 
(replacing Ti2 in the subsurface layer) are shown, respec- 
tively. 

For the first case of adsorbates shown in Fig. 4(a), the 
energy gap is full of six defect states, five of which are 
occupied and one that is close to the conduction-band 
is empty. Fermi level is at 1.84 eV. The system has a 
narrow gap of only 0.15 eV. The gap is indirect and the 
highest occupied defect state makes a maximum at mid- 
way, otherwise almost flat-band, along XM. The lowest 
unoccupied band is also a defect state making a minimum 
at r and having the shape of the lowest conduction-band 
along most directions except around A'. Its LDOS has a 
two-dimensional character as expected, and is a result of 
antibonding interactions between electrons from Pt and 
01 ions. The lower group of four bands all have Pt con- 
tributions with three upper bands having contributions 
from 01. Especially, the flat-band at 1.0 eV causing a 
peak in the LDOS and the peak above that are mainly 
due to 01 ion, whereas the lowest of the four has con- 
tributions mainly from 02 ion. This band and the one 
below, separated by a pseudo-gap, which is barely above 
the valence-band top, are of the same character, namely, 
of Pt and 02. 

The second adsorbate case described in Fig. 3(b) is en- 
ergetically more preferable than the first one as they both 
contain equal number of atomic species. In this case [see 
Fig. 4(b)] the same five defect states fill the bulk band 
gap as in the previous case. All of them lie below the 
Fermi level which is at 1.95 eV relative to valence-band 
top and the gap is 0.13 eV, again small. The lowest unoc- 
cupied state being bulk-like this time is the conduction- 
band minimum at T. The empty defect state is higher 
in energy. The highest occupied defect state has a max- 
imum at X' , however the band is almost flat along AT 
making the gap nearly a direct one. The band is also flat 
in some part of the way along TM, both causing higher 
densities in the LDOS picture. The contribution by the 
01 ion is mainly in the upper defect bands in the energy 
range between E F and about 0.2 eV below that along 
with Pt ion, while that by 02 is mainly in the lower de- 
fect states starting from the highest peak at about 0.6 
eV down into the valence-bands. 

In the case of surface substitutional, the bulk energy 
gap is full of four defect states as seen in Fig. 4(sl), 
lower two of which are fully occupied and the upper two 
are almost half-filled each. Therefore the system shows 
metallic behavior. Fermi level is at 1.51 eV relative to 
the valence-band top. These two upper bands cross each 
other along high-symmetry directions of the surface Bril- 
louin zone at three places, namely, very close to the T 
point along AT, and passed midway along the TM and 
XM directions. The LDOS due to Pt ion as well as some 
due to neighboring two oxygen ions 01 and 02 are com- 
ing from these bands dispersed around the Fermi level. 
The two lower bands are surface states due to oxygen 
ions 01 and 02. 
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Figure 4(s2) shows the band structure for the subsur- 
face substitutional where the energy gap is full of four 
defect states that are all occupied. Fermi level is at 1.85 
eV relative to the valence-band top giving rise to an in- 
direct energy gap of 0.23 eV. The LDOS of the highest 
occupied defect band having a sharp peak due to the al- 
most flat part along XT is contributed by Pt ion and 
its neighboring oxygen ions. This band is pretty similar 
in shape to the half-filled one in (si) case that makes 
a maximum at M. Next piece of LDOS between 0.3 
and 0.9 eV is due to the surface oxygen 01. This band 
was observed to be similar to the clean surface band due 
to 01 in Fig. [2fa), especially around M point along all 
three directions. The band that is hardly visible above 
the valence-band top at M point is also similar to the 
second band of clean surface at the same region. Since 
the Pt ion is deeper than skin of the surface region in 
this case the surface is heeled and the surface states are 
back. The van Hove singularity at about 0.9 eV, where 
this band ends at the right-most M point in Fig. 0|b), is 
also evident just like the one at Ep=QA0 eV in Fig. [2a). 



B. Dilute impurities 

TiO2(001)-(2x2) unit cell prevents charge transfers be- 
tween the impurity sites by providing a 7.53 A Pt-Pt sep- 
aration together with local screening effects through the 
slab. Therefore, both the geometric and the electronic 
structures show significant differences from the (lxl) 
counterparts. For instance, Pt interstitials can be stable 
inside (2x2) slab without segregation effects as opposed 
to the case with (lxl) surface. For the incorporation of 
Pt ion as an adsorbate or as an interstitial, we examined 
a total of six cases that are shown in Figs. EJa)-[5jf ) , in 
the first three of which the Pt is situated on/in the sur- 
face and in the last three the Pt ion is below the surface 
layer. 

In the first case, [Fig. 03a)], Pt adsorbate is twofold 
coordinated with nearest-neighbor Ti5c and 02c along 
the Ol-Til row. Pt raises the 02c upward by 0.92 A and 
pushes it in [100] direction by 0.83 A from its relaxed 
lattice position. As a result of the SMSI, bond distances 
from Ti5c to Pt, Pt to 02c, and 02c to the next Ti5c, 
which are connected successively through line segments 
in a row, become 2.36, 1.90, and 180 A, respectively. The 
other Ol-Til row remains to be less effected by the pres- 
ence of Pt adsorbate with Ti5c-02c bonds, 2.11 and 1.81 
A in length, making a Ti5c-02c-Ti5c angle of 144.4° as in 
the case of relaxed clean surface. Among the side bonds, 
other than the skewing toward [100], merely the Til-03 
one, which connects to the promoted 02c, is elongated 
by 0.04 A with a dihedral angle of 7.23° bearing a subtle 
difference from the values calculated for the (2x2) clean 
surface. Moreover, the distances and angles are not dis- 
torted considerably in the second Ti02 layer near the 
surface. 

The structure shown in Fig. ED(b), being the lowest 



energy configuration among adsorptional and interstitial 
cases, undergoes a more complex atomic rearrangement 
upon Pt adsorption. The adsorbate migrates to a bridge 
position in between Til and 02 with bond lengths of 
2.78 and 2.07 A for Pt-Til and Pt-02, respectively. As 
a consequence of the SMSI, the nearest-neighbor 02c is 
pushed from its lattice site upward above the midpoint 
between Til and Pt forming a triangle. It is aligned 
perpendicular to the surface, with 02c-Til and 02c-Pt 
sides of 1.81 and 1.95 A, respectively. In this geometry 
02c elevation above the Pt adsorbate is calculated to be 
0.83 A. Another triangle forms between the consecutive 
Til on the Til-Ol row, neighboring 02 and Pt which 
lies at the vertex connecting the two triangular atomic 
arrangements in a fourfold coordination with these sur- 
face atoms. In this second triangle Pt-Til bond distance 
becomes 2.30 A. Furthermore, the two Ti5c-02c bonds 
become slightly distorted, 1.81 and 2.17 A in length, with 
Ti5c-02c-Ti5c angle of 143.9° on the second Til-Ol row. 
The only noticeable difference occurs in the Til-03 bond 
which extends to 2.16 A. All other side bonds preserve 
the skewing similar to the (2x2) clean surface results. 

The third adsorptional model adopts a more neatly 
symmetrical structure among Pt/TiO2(001)-(2x2) cases 
as shown in Fig. Etc). Pt resides at the midpoint between 
two undercoordinated surface oxygen with a bond length 
of 1.99 A. Meanwhile, Pt also makes equidistant bonds 
with four Ti5c's, each of which 2.8 A in length. This 
extended bond distance (see Table [TTJ) entails a rather 
weak Pt-Ti5c interaction on the surface. As a result of 
this isotropic fourfold coordination with surface Ti's, the 
axial Til-03 bonds become aligned parallel to [001] direc- 
tion. On the other hand, the two Ti2-02 bonds, being 
coplanar with the Pt, are the only slanted ones which 
make two dihedral congruent angles of 12.3° with the 
(100) and (100) planes. 

When Pt is substituted for one of the surface Ti ion 
(Til), the Pt-03 bond, although unaltered in length, no- 
ticeably differs from Ti-03 bonds by skewing more from 
2.21° to 18.21°. Due to the excess charge incorporated by 
the substitutional impurity, the Pt-02-Til row exhibits 
slight distortions in bonds and angles. Resultant two Pt- 
02c-Ti5c angles become 141.8/153.5° with Pt-02c and 
Ti5c-02c bond distances of 1.92/2.10 and 2.01/1.82 A, 
respectively. On the other hand, the Til-02-Til row re- 
tains the values as those of the clean surface. 

As in the case of (lxl), Pt substitution for Ti2 re- 
sults in the reduction in the skewness in the axial bonds 
as they align parallel to [001] direction. While the axial 
bond lengths assume an ideal-like value of 2.02 A, the 
nearest-neighbor equatorial Pt-0 distances become 1.99 
A which describes a slightly extended value over the Ti- 
O equatorial bond length of 1.94 A. Therefore, the only 
substantial displacement is obtained for the atoms at the 
surface Ti02 layer along [100]. This rearrangement re- 
duces the skewness but preserves the non-symmetric na- 
ture of Ti5c-02c bonds which read 2.14 A and 1.81 A 
with a Ti5c-02c-Ti5c angle of 144.8°. In other words, 
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TABLE II: Calculated values of some key parameters for Pt-TiO2(001) anatase system: work function, Fermi energy relative 
to bulk valence top, and change in energy-band gap, as well as Pt-depth relative to surface oxygen, Pt-O distance, and Pt-Ti 
distance for each model. Labeling of models follow from Figs [3] and [5] for (lxl) and (2x2) surfaces, respectively. 



Surface 


Model 


W(eV) 


E F (eV) 


A£g(eV) 


hpt(A) 


dpt-o(A) 


dpt-Ti(A) 


lxl 


Clean 


6.88 


0.40 


0.00 










(si) 


5.83 


1.51 


_ 


-0.46 


1.94(02c),1.96(02),2.18(03) 


_ 




(s2) 


5.36 


1.85 


-1.45 


-3.04 


1.95(O3,O4),2.02(O2,O5) 


_ 




(a) 


4.16 


1.84 


-1.53 


0.76 


1.97(01) 


2.55(Til) 




00 


5.21 


1.95 


-1.55 


0.54 


1.96(01) 


2.61(Til), 2.71(Til) 


2x2 


Clean 


6.77 


0.38 


0.00 










(si) 


6.13 


1.12 


-1.13 


-0.68 


1.92(02c), 2.10(O2c), 1.99(03) 


— 




(s2) 


6.25 


0.88 


-0.50 


-3.04 


1.99 (03,04), 2.02 (02,05) 






(a) 


5.86 


0.96 


-0.58 


1.24 


1.90(01) 


2.36(Til) 




(b) 


5.61 


1.41 


-1.03 


0.24 


1.95(Ol),2.07(O3) 


2.30(Til),2.78(Til) 




(c) 


5.64 


1.80 


-1.42 


0.02 


1.99(01) 


2.80(Til) 




(d) 


5.71 


1.51 


-1.13 


-2.91 


l,97(O3),2.06(O3) 


2.72(Til),2.63(Ti2),2.86(Ti2) 




(e) 


5.69 


1.49 


-1.11 


-3.59 


2.01(04) 


2.71(Ti2),2.77(Ti2),2.91(Ti3) 




(f) 


5.62 


1.56 


-1.18 


-5.29 


2.01(05) 


2.88(Ti2),2.75(Ti3) 



as the number of Pt substitutions for four possible Ti2 
ions increases, the structure adopts an ideal-like geome- 
try which gains a plane mirror symmetry as in the case 
of s2-lxl. 

For anatase TiO2(001)-(2x2) surface, being contrary 
to the (lxl) cases, Pt is found to be stable once it pen- 
etrates into the interstitial cavities, starting from the 03 
level, as shown in Figs. [5^d)-[5jf ) . Pt ions are encapsu- 
lated by the slab in an octahedral position at the mid- 
point between the two level oxygen in interaction with 
the six nearest-neighbor Ti's. The location of the en- 
capsulated Pt is determined dominantly by the strong 
Pt-0 interaction, as being at the same depth with and 
in the middle of, the fully coordinated two level oxygen. 
This leads to two structural ramifications as fingerprints. 
First, two oxygen atoms whose interconnecting line is 
perpendicular to O-Pt-0 bonding, in the closest (pre- 
ceding or succeeding) oxygen atomic layer, are slightly 
repelled out from their lattice positions as a result of the 
induced stress due to the excess charge brought by the 
impurity. This is not solely specific to interstitial cases. 
We obtained the very similar result for the adsorptional 
model in Fig. [5jc). Second, Pt ion, making equidistant 
equatorial bonds with, maintains four Ti's at the corners 
of a square shape which can be seen through [001] di- 
rection. This is due mainly to Pt impurity rather than 
the TiC>2 lattice itself since the third adsorptional case 
[see top view of Fig. [U^c)] adopts surface Ti's in the same 
geometry. This is so that Til-03 axial bonds are aligned 
parallel to [001] as opposed to their slanted posture in the 
clean surface. Evidently, similar arguments apply for the 
first two interstitial cases [Figs. [5^d) and[5je)], as well. 
This time Pt holds Ti2's in a squarely manner causing 
Ti2-02 bonds to align vertically with respect to the sur- 
face plane. When Pt penetrates deeper than the second 



Ti02 layer into the cavities, the symmetry breaking over 
the Ti5c-02c bonds and the skewness of Til-03-Ti2-02 
ring, similar to the clean surface can finally be repro- 
duced. This corresponds to the last interstitial model as 
shown in Fig.^f). Therefore, Pt implants induce a local 
stress causing nearby atoms to slightly rearrange from 
their lattice positions in the TiC>2(001)-(2x2) surface as 
a result of the strong Pt-0 interaction. 

The electronic structure for the geometries in 
Figs. Efa)-r5][f) with the interstitial impurity are pre- 
sented in Figs. [Ha)-[H]Jf), respectively, and moreover in 
Figs. Efsl) and[B]Js2), those for the substitutional impu- 
rity at the surface (replacing Til) and that in the bulk 
(replacing Ti2 in the subsurface layer) are shown, respec- 
tively. 

For the first adsorbate case described in Fig. Ufa), the 
impurity bands due to the interaction between Pt and 
the neighboring promoted-oxygen 01' (prime for bonded 
neighbor) are grouped in three sets within the bulk band- 
gap region [see Fig. E^a)]. Two of them being almost 
conjugate of each other are in the energy ranges of 2.3- 
2.6 and 0.15-0.5 eV, respectively, while the bands in the 
third group, being almost flat, producing a sharp peak 
in the LDOS picture, are located in between the other 
two groups at around 1 eV. This is the highest occu- 
pied set of states of mainly Pt-Ol' character with some 
contributions from the 02s. The Fermi level is at 0.96 
eV, and since the lowest unoccupied band is the bulk 
conduction-band having a minimum at T the energy gap 
is direct and 1.12 eV in width. The conjugate bands 
are dispersed along k y but fiat along TJ and KJ'. The 
empty one is in resonance with the bulk bands around 
r J and localized in the gap otherwise. Due to this asym- 
metry in localization the LDOS is also not symmetric in 
shape. The filled one being in the gap has a symmetrical 
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FIG. 5: (Color on-line) Interactions of single Pt atom with anatase TiC>2(001) surface (in order of penetration depth of the Pt 
atom): (a) Pt adsorbed on the surface, (b) surface oxygen supported by Pt with SMSI, and [(c)-(f)] stable structures of Pt in 
the surface. Pt atom is shown in white and denoted by the biggest sphere while Ti is in gray and O is represented as small red 
(dark) balls. 



one-dimensional (ID) LDOS coming from Pt-Ol' inter- 
action, with a smaller mixture of Pt and larger mixture 
of 01', which is the other way around in the empty con- 
jugate state. In addition to bonding to 01', Pt causes 
slightrcpositioning to other surface oxygen, as a result 
of which, a little contribution comes from 01 along the 
same line as Pt-Ol' [see Fig. 4(a)] in the upper part 
of the delocalized ID-peak. The lower peak of the same 
LDOS is degenerate with a contribution due to a new flat 
band of 02 character adjacent to the edge of the above 
band. The lower two bands are surface-like and again of 
02 characters. 

The case shown in Fig. [5](b) is the most favorable one 
with the lowest total energy among all (except the substi- 
tutional cases since they belong to different stoichiometry 
which prevents a direct comparison of their total ener- 
gies) and its electronic nature is again a semiconductor, 
like all the other dilute impurity cases. Comparing its en- 
ergy bands [see Fig. H^b)] with the case of Fig. [^a), the 



empty impurity band (not shown) is now pushed up into 
the bulk continuum of the conduction-band. The two flat 
peaks (Pt-Ol'), split by 0.15 cV, arc also raised in energy 
causing the Fermi level to be at 1.41 eV, and decreasing 
the (direct) energy gap to 0.67 eV. They have also contri- 
butions from 01 and 02'. The next lower band (due to 
Pt-01'-Til-02 chain) is now much less dispersed and it is 
separated from the impurity band below, being no longer 
degenerate at T as in case of Fig. EJa). The fifth (shown) 
band is also due to Pt-Ol' and Til-Ol interactions. 

Figure[5Jc) shows the last of the adsorbate cases where 
Pt is situated very symmetrically on the 01 layer. All 
four Til's are equivalent; two of the four Ol's are bonded 
to the Pt interstitial, and two of the four 02 's are dis- 
placed as seen in top view of Fig. [5jc). Consequently, 
there are eight bands fallen into the gap. The first one 
from the top, again empty, is within the conduction- 
band, with localization around K, and conjugate to the 
sixth band. It is mainly of Pt character. The second 



10 




(a) 













r 

(b) 













































f 




• 


























































t 





FIG. 6: (Color on-line) (si) Energy bands for the anatase 
TiO2(001)-2x2 surface with Pt adsorbed as a substitutional 
impurity at a surface-Ti site, (s2) at a subsurface- Ti site, and 
[(a)-(f)] for the interstitial impurity cases in Figs. [5Ja)-[SJf). 



band is the highest occupied Ol'-related surface band 
which makes a maximum at T and flat along TJ. The 
Fermi level is at 1.80 eV and the (direct) gap is 0.28 eV, 
the narrowest gap of all cases. The third band at about 
1.5 eV is an impurity band due to Pt-02' interaction of 
larger Pt mixing, whereas the fourth band, at about 1 
eV, has a sharp peak of same character with larger 02' 



mixing. The fifth band is a perfect flat one due to dis- 
placed 02s since Pt pushes 02s off their lattice positions 
causing these stress-induced flat going states that reflect 
major oxygen character due to relatively weaker neigh- 
boring interactions. The sixth band is a rather dispersed 
two-dimensional surface band of 01', 02, and 01 charac- 
ters, in decreasing order of contribution, which is almost 
the mirror image of the empty first band above. They 
are both very symmetric with respect to k x and k y as 
expected [see top view of Fig. [S^c)]. The last two, the 
seventh and eighth bands are impurity bands of mainly 
02', 02, and slightly of 01' and Pt contributions. 

Since the geometric structures, shown in Figs. [5^d)- 
[5jf), of the Pt ions as bulk interstitials (subsurface and 
deeper) in their nearest-neighbor environment are equiv- 
alent, their electronic band structures are also very sim- 
ilar, especially for (d) and (f). It is expected that as the 
Pt ion is placed deeper into the bulk [see Figs.|5(d)-|Bff)], 
the band structure will alternate as "(d) and (e)". The 
empty defect band lies partly in the gap. The second im- 
purity band being filled is due to Pt-01' (02, 03', 04') 
interactions, and making a maximum at J, J', and J in 
cases (d), (e), and (f), respectively. Similarly, the Fermi 
level is situated at about 1.5 eV and the energy gap is 
about 0.5 eV (see Table II). The lower filled bands are 
the mixture of surface bands and less dispersed impurity 
bands due to Pt interstitials. 

When the surface Ti is substituted by a Pt impurity on 
a 2x2 reconstruction, the gap is filled by several impu- 
rity bands as seen in Fig.^sl). One of them is far above 
in the conduction-band region and partly fallen into the 
gap. The next lower one is an empty band making a 
minimum at J' point around ~1.6 eV. The bandwidth 
is about 0.2 eV due to the interaction between the Pt 
substitutional impurity with its neighbors 03' and 01'. 
Below is the highest occupied band of surface-like state 
dispersed by the interaction between surface Ti and 02'. 
This band makes a maximum, Ep — l.Vl eV, at J' point 
as well, causing the 0.57 eV gap to be a direct one and 
narrowed as compared to the 1.70 eV gap for the clean 
surface. The band is nearly flat along JK and K J' which 
results in a sharp peak in the LDOS at Fermi level. Hav- 
ing the upper empty band going almost flat along KJ', 
as well causing parallel bands in this direction, one may 
expect an enhancement in the optical transitions rate. 
The two bands further below are due to 03'-Pt-01' and 
Ti-02' bonding interactions, respectively. Morover the 
bands below those look more like, roughly, the surface 
states of clean (001). 

In the case of placing the substitutional Pt ion at the 
second layer we have a bulk-like impurity problem where 
the surface layer is chemically similar to the clean surface 
with four Ti ions exposed. This is also evident in the 
band structure shown in Fig. [S£s2). The empty states 
fallen from the conduction-band into the forbidden gap, 
but still above the conduction-band minimum, are very 
symmetric with respect to k x and k y directions. And in 
the optical gap, most of the bands look surface-like with 
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impurity bands passing through them. The Fermi level is 
at 0.88 eV and the gap of 1.20 eV is again a direct one at 
r like in the clean surface case. The difference of 0.50 eV 
is partly (about 0.30 eV) due to the un-reconstruction 
of the clean surface with the subsurface impurity sub- 
stitution. The Pt based impurity state is almost flat at 
around 0.2 eV above the valence-band top. 

The strong dispersion of defect bands seen in (lxl) 
cases and the rather less dispersed nature of these bands 
in (2x2) models originating from the Pt impurities sig- 
nify the role of Pt-Pt interaction in relation to Pt concen- 
tration. Furthermore, the distinguishable flatness of the 
impurity bands shown in Fig. (6^a) and ^h) arise from 
the minimal coordination of Pt adsorbate with the sur- 
face ions as a result of its spatial location, in addition to 
the diluteness of these impurities. 

1 ML Pt substitution for undercoordination drive 
anatase into a metallic state. At the same coverage Pt 
adsorption, and substitution for the second-layer Ti in 
this surface yields a low band-gap semiconducting sys- 
tem which would be active in the infrared region. On 
the other hand, anatase TiO2(001) can be functional- 
ized for visible activity by Pt impurities implanted with 
a (2x2) periodicity for all models except the (c) case 
[see Fig. EJc)] which corresponds to a band gap narrow- 
ing of 1.42 eV. For this impurity concentration, this rep- 
resents the maximum value close to those of the 1 ML 
cases. LDOS analysis indicates that the number of de- 
fect states derived from the valence bands increase with 
the coordination of Pt with 01 and 02 ions which is 
also maximum due to its high-symmetry relaxed posi- 
tion. Therefore, even though, platinized Ti02 is known 
to give visible-light activity^ our calculations show that 
band gap narrowing depends on the impurity concentra- 
tion and the coordination number of Pt with the near 
surface oxygen. 



C. Analysis of the electronic density 

The topological analysis of the charge density gives ac- 
curate information about the bonding characteristics for 
the neighboring atoms. Therefore, one needs a qualita- 
tive description of the interatomic charge distributions 
which can be computed by employing Bader analysis 
based on atom in molecule (AIM) theory. To do so, the 
real-space cell is partitioned into Bader volumes delim- 
ited by local zero-flux surfaces of the electron density 
gradient vector field. Then, these volumes can be in- 
tegrated around an atomic region to calculate the local 
charge depletion and accumulation. We obtained these 
atomic properties using the AIM formalism that is imple- 
mented with a recent grid-based algorithm4£ The Bader 
charge results for Pt and surface region Ti and O atoms 
are presented in Table [TTTl in the cases of (1 x 1) and (2x2) 
unit cells. For (2x2) surface, there are four possible Ti 
(or O) ions to choose from which are at the same atomic 
layer. We have preferred to provide the calculated values 



TABLE III: The valence charge accumulation based on Bader 
analysis for Pt-TiO2(001) anatase system. Atom labels follow 
Fig. [1] The charges for lattice atoms that are closest to Pt 
are presented in the case of (2x2) slab. 



Surface 


Model 


Ol 


02 


03 


Til 


Ti2 


Pt 


lxl 


Clean 


-1.26 


-1.34 


-1.33 


+2.61 


+2.64 


— 




(si) 


-0.79 


-0.99 


-1.15 


— 


+2.62 


+1.63 




(s2) 


-1.27 


-1.20 


-1.02 


+2.57 


— 


+1.84 




(a) 


-1.15 


-1.34 


-1.32 


+2.52 


+2.64 


-0.03 




(b) 


-1.17 


-1.33 


-1.31 


+2.39 


+2.64 


+0.11 


2x2 


Clean 


-1.26 


-1.34 


-1.33 


+2.60 


+2.64 


— 




(si) 


-1.07 


-1.07 


-1.33 


+2.58 


+2.63 


+1.64 




(s2) 


-1.26 


-1.19 


-1.18 


+2.60 


+2.64 


+1.85 




(a) 


-1.07 


-1.35 


-1.33 


+2.49 


+2.64 


+0.02 




(b) 


-1.08 


-1.31 


-1.33 


+2.57 


+2.65 


+0.06 




(c) 


-1.22 


-1.32 


-1.33 


+2.54 


+2.64 


+0.07 




(d) 


-1.26 


-1.33 


-1.23 


+2.56 


+2.56 


+0.08 




(e) 


-1.23 


-1.34 


-1.27 


+2.58 


+2.55 


+0.15 




(f) 


-1.26 


-1.34 


-1.31 


+2.60 


+2.59 


+0.14 



for the ones that are closer to the Pt impurity, in order to 
make a better comparison with the (lxl) counterparts. 




FIG. 7: (Color on-line) 3D charge-density plot for the "b- 
2x2" case of Pt on anatase TiO 2 (001) surface [Fig. Ob)]. 

In the case of clean surface, a deep lying Ti ion, which 
should reflect bulk-like properties, transfers an amount of 
0.443 electronic charges to each of the neighboring oxygen 
losing its last atomic shell. Hence, this fully coordinated 
Ti ion accumulates a Bader charge of +2.66e. A bulk-like 
oxygen, on the other hand, gets a valence charge state of 
— 1.33e. Calatayud et aL— computed these values for the 
bulk anatase Ti02 as Qxi = 2.96e and Qo = — 1.48e us- 
ing a different exchange-correlation scheme. The Bader 
charge accumulation around the ions is also sensitive to 
the determination of integration regions with boundaries 
delimited by zero gradient of the electronic density. Al- 



12 



though we obtained slightly lower charge states for Ti 
an O ions resulting in a relatively less polarized bond- 
ing, in quite an agreement both results obey the same 
stoichiometry by correctly adding up to charge neutral- 
ity of TiC>2 and are smaller than the nominal oxidation 
states obtained for a generic ionic oxide such as MgO. As 
a result, a polarized covalent bonding develops between 
charged Ti cations and O anions. 

Naturally, the AIM charge values, presented in Ta- 
ble HIH for (lxl) and (2x2) periodicities in the case of 
clean surface models exhibit a one-to-one correspondence 
as expected, since the electronic properties (such as the 
band gap, work function, etc.) derived from their charge 
densities must represent the same surface. The oxidation 
states of surface layer ions affect the reactivity of single 
crystals of anatase (001). Hence, we computed the Bader 
charges for the undercoordinated ions as — 1.26e for 01 
and +2.61e for Til being lower by ~ 5% and ~ 1% than 
that of the deep lying O and Ti, respectively. More- 
over, bulk-like charge states are adopted starting from 
02 atomic layer which stays 2.56 A below the surface 
oxygen. Clearly, bulk termination bears slight differences 
in charge states of deep and surface ions, particularly in 
the case of surface oxygen, indicating a rather low reac- 
tivity of the clean (001) surface. This prediction is in 
good agreement with the experimental observations that 
clean surfaces of Ti02 exhibit lower catalytic activity 
than stepped (101) and oxygen-defect (001) surfaces^ 

Pt incorporation yields significant disturbance in the 
electronic density in the vicinity of the impurity site. Par- 
ticularly, this effect is observed for the oxygen in the close 
proximity of Pt while the electron depletion from around 
the nearest-neighbor Ti ion remains minimal upon Pt de- 
position. From the data presented in Table [TTTl the stan- 
dard deviation in the Bader charges with respect to the 
reference clean surface values have been calculated to be 
0.17, 0.14, and 0.12 for 01, 02, and 03 while the corre- 
sponding values arc 0.09 and 0.04 for Til and Ti2, respec- 
tively. Smaller deviations obtained for the deeper lying 
ions also imply a limited contribution of these atoms to 
the DOS and surface bands near the Fermi level. There- 
fore, for instance, we obtain the smallest band gap for the 
adsorptional case (c) among (2x2) structures although 
it has lower number of Pt-0 interactions than interstitial 
configurations do. 

When substituted for Ti ion either at the surface 
or in the slab, Pt can not acquire the same oxidation 
state that Ti had. This results in a relatively lower 
amount of charge accumulation around the neighboring 
oxygen resulting in a less polarized covalcncy between 
Pt and O. This is clear also, for instance, in the three- 
dimensional (3D) charge-density plot presented in Fig. [71 
which belongs to the case (b) of Pt-Ti02 (2x2) system. 
Ti's show lesser valence electronic density distributions 
around them, indicating strongly polarized covalent Ti- 
O bonding. The electron depiction from Pt to 01 is 
lower than, that, for instance, from Ti2 to neighboring 
oxygen. Moreover, the charge density around 01 is no- 



ticeably smaller than that of 03 (and also of 02). In 
summary, Ti-0 bond polarization, and therefore interac- 
tion, in Ti02 lattice environment is stronger than that 
of Pt-O. 

On the other hand, the strength of Pt-Ti02 interaction 
can also be compared with respect to different model 
cases presented in Table Mil Bader charges calculated 
for Pt ion suggest that it interacts with the lattice in the 
substitutional cases stronger than in the adsorptional and 
interstitial configurations. 

D. Thermodynamic stability of the phases 

The (lxl) and (2x2) surface supercells comprise un- 
equal amounts of atomic species. Moreover, adsorp- 
tional or interstitial cases, being stoichiometrically differ- 
ent from the substitutional ones, represent an addition of 
impurity instead of a Pt-Ti replacement. When the su- 
percell total energies are considered, the lowest energy 
structures are b-lxl [Fig. [j^b)] and b-2x2 [Fig. GJb)] 
among the adsorptional/intcrstitial cases while they turn 
out to be sl-lxl [Fig. EJsl)] and s2-2x2 [Fig. [5(s2)] for 
the substitutional geometries at one and quarter ML con- 
centrations, respectively. Therefore one cannot directly 
compare their stability by just looking at their supercell 
total energies. 

We employed formalism of Qian et a/42, and 
Northrup4^ to study the thermodynamic stability of the 
Pt-incorporated TiO2(001) surfaces that have varying 
number of constituents at different concentrations. In 
this scheme, relative formation energy is defined as a 
function of the chemical potential of the excess atomic 
species as 

Eform = -Ept/Ti0 2 ~ ^TiOs. ~ Anxi^Ti — Anpt/Upt 

where -Ept/Ti0 2 an d ^Ti0 2 are t ne total energies of the Pt 
incorporated and bare Ti02 surfaces while /ixi and /xpt 
stands for the chemical potentials of Ti and Pt. An-pi and 
Anpt represent the differences in the number of atoms of 
each atomic species with respect to the reference clean 
surface. Formation energy in this form is a function of 
Ti and Pt chemical potentials. Equilibrium is reached 
when the chemical potential of a given species is equal in 
all the phases that are in contact with each other. Also, 
these phases must be in equilibrium with bulk anatase 
such that 

MTi + 2/io = MTi0 2 • 

This relation interrelates chemical potential of Ti to O 
chemical potential that varies accordingly with the ex- 
perimental conditions. The value of (iTi must be smaller 
than that of the hep Ti bulk solid phase which is an 
undesirable formation at the surface and is referred as 
the Ti-rich conditions. The other extreme is that when 
the surface oxygen are found in thermodynamic equilib- 
rium with the molecular oxygen bath corresponding to 
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(a) Thermodynamic stability of Pt/TiO 2 (001) surfaces 
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(b) Phase diagram of Pt/TiO 2 (001) surfaces 



I -I 



— 




O-rich 



Ti-rich 



FIG. 8: (Color on-line) (a) Normalized formation energies 
of Pt-incorporated TiO2(001) structures relative to that of 
clean surface as a function of the Ti chemical potential. (For 
numerical figures, see Table HV[) Pt impurities are chosen to 
be in thermodynamic equilibrium with fee Pt bulk phase, (b) 
Phase diagram of Pt/TiO2(001) surface as a function of Pt 
and Ti chemical potentials. 



O-rich conditions. Assuming it as an ideal gas, the in- 
termolecular interactions can be neglected. Hence, the 
chemical potential of oxygen can be referenced to the 
value at the O2 molecule which is no = Eo 2 /2. There- 
fore, Ti chemical potential relative to its bulk value 
A/ixi = MTi — MTi var ies between A/iTi = (Ti-rich 
conditions) and = —8.9 cV (O-rich conditions), for 

Ti02- Since the relative formation energy of phases is 
defined as a function of Ti and Pt chemical potentials, 
Pt is assumed to be in thermodynamic equilibrium with 
its fee bulk solid phase that serves as a reservoir for Pt 
atoms. Hence, [j-Pt is chosen to represent such an ex- 
perimental condition. In fact, this corresponds to just 
another extreme for undesircd Pt phases on the surface 
that can be avoided by [ipt < Mpt' k - 

The formation energies are shown in Fig. [8^a) relative 
to that of the clean surface for 12 phases as a function 



TABLE IV: Numerical values for the formation energy of 
Pt/Ti02 surfaces for adsorptional and interstitial phases (all 
in eV/lxl) referring to the constant valued functions de- 
picted in Fig. [5] Relative formation energies are given with 
respect to that of the bare surface. Those of the substitu- 
tional phases not provided since they vary with varying Ti 
chemical potential due to different surface stoichiometry. 



Phase 


-E-form (Absolute) 


Sform (Relative) 


Clean 


0.76 


0.00 


a-lx 1 


2.76 


2.00 


b-lxl 


2.53 


1.77 


a-2x2 


0.59 


-0.17 


b-2x2 


0.39 


-0.37 


c-2x2 


0.53 


-0.23 


d-2x2 


0.94 


0.18 


e-2x2 


0.87 


0.11 


f-2x2 


0.87 


0.11 



of Ti chemical potential over the full range of its allowed 
values. They are normalized to a lxl unit cell in order 
to compare their thermodynamic stabilities. Pt adsorp- 
tion at 1 ML coverage results in unstable surfaces with 
relative formation energies at 1.77 and 2.00 eV/lxl for 
(b) and (a) adsorption modes, respectively. Substitu- 
tional cases at this concentration are the most unstable 
phases within the range of —8.56 < A/ixi < eV with in- 
creasing instability toward Ti-rich conditions while they 
become the most stable phases under O-rich conditions 
for —8.90 < A/iTi < —8.56 eV in favor of substitution 
for Ti5c (sl-lxl). For this latter experimental situation, 
in which the surface is in thermodynamic equilibrium 
with the molecular oxygen, s2-lxl phase has a relatively 
higher formation energy by 0.10 eV/lxl followed by s2- 
2x2 and sl-2x2 structures which are energetically less 
favorable by 0.26 and 0.33 eV/lxl. Similarly, over the 
range -7.15 < Afj, Ti < eV, Pt substituted TiO 2 (001)- 
(2x2) surfaces are the most unstable structures among 
0.25 ML phases whereas they turn out to be the most sta- 
ble cases, in which the formation energy is slightly lower 
for s2-2x2 phase, over the range —8.56 < A/xxi < —7.15 
cV closer to O-rich conditions. 

On the other hand, Pt interstitials are unstable rel- 
ative to the formation energy of the clean surface by 
0.18 eV/lxl for (d), and by 0.11 eV/lxl for (e) and 
(f) at (2x2) reconstruction implying increased stability 
with increasing Pt penetration depth. These results sug- 
gest that practical interstitial applications might require 
thermal treatment. Therefore, our calculations for the 
stabilities of substitutional cases under O-rich conditions 
and of Pt interstitial phases are in good agreement with 
the experimental results of Zhang et al£- who suggested 
that neutral Pt atoms can thermally diffuse into Ti02 
lattice under oxidizing atmosphere. They also argue that 
these diffused Pt atoms can either substitute for the Ti 4+ 
sites when oxidized to Pt 2+ (for which our calculations 
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show a charge state of +1.85e) or they occupy interstitial 
sites. 

In contrary to those of the 1 ML phases, Pt adsorbates 
at 0.25 ML happens to be energetically more stable rel- 
ative to bare surface by 0.17, 0.23, and 0.37 eV/lxl for 
the (a), (c), and (b) adsorption phases, respectively, for 
the whole range of allowed Ti chemical potential. Be- 
sides, b-(2x2) structure is also the most stable phase 
within —7.15 < A/iTi < eV which spans 80.3% of the 
whole range from Ti-low to Ti-rich conditions. 

The phase diagram shown in Fig.[8jb) has been derived 
from the results obtained for the energetically more sta- 
ble phases, presented in Fig. [5{a), for varying Ti and Pt 
chemical-potential values. Hence, high formation energy 
surfaces have not been considered due to their instability. 
Under Pt-rich conditions, which refer to the formation 
energies presented in Fig.OJa), three most stable phases 
exist for varying chemical potential of Ti. Under Pt-rich 
and O-rich conditions the most stable phase is sl-lxl in 
which all of the surface-Ti ions substituted with Pt. A 
small deviation from these conditions by slightly increas- 
ing the Ti chemical potential switches the phase to 0.25 
ML concentration surface of s2-l x 1. For lower O concen- 
trations b-2 x 2 surface is more stable and dominant over 
the range of allowed Ti chemical potential. Under Pt- 
poor conditions the phase diagram reproduces the clean 
TiC>2(001) surface with no impurities. 



IV. CONCLUSIONS 

We systematically studied the structural and electronic 
properties of Pt impurities in the form of adsorptional, in- 
terstitial, and substitutional cases for anatase TiC>2(001) 
with (lxl) and (2x2) surface periodicities. The former 
represents full coverage while the latter corresponds to 
isolated impurities. Depending on the Pt concentration 
per unit cell area, impurity-impurity electron coupling 
strength mediates the mode of atomic rearrangements as 
they are clearly different for (lxl) and (2x2) models. 
For instance, Pt adsorption at the bridge site obtained 
for 1 ML coverage in Fig. [3{a) is corresponded by the 
pattern shown in Fig.[5{a) at 0.25 ML coverage, in which 
surface oxygen is promoted by the adsorbate as a result 
of the SMSI. This difference is even more pronounced 
for interstitial cases. When implanted inside the slab for 
full coverage, Pt atoms form parallel metallic wires inside 
TiC>2 where interlayer distances slightly increase due to 
local segregation, while Pt impurities can be encapsu- 
lated by the (2x2) lattice at interstitial cavities to form 
structures without undergoing a major reconstruction. 

In addition, another dominant factor in the formation 
of low energy Pt/TiC>2 structures is the nearest- neighbor 
Pt-0 coordination which derives from the impurity- 
lattice oxygen charge transfer. Relative stabilities of 
these structures can be addressed to local disturbance 
on the potential-energy surface induced by the excess 
electrons brought by the impurities that consequently ac- 



count for the enhancement of the electron trapping effi- 
ciencies. 

Clean TiO2(001) possesses surface states derived from 
the valence-bands in the energy-band gap near the VBM 
originating from the undercoordinated surface oxygen. 
The nature of these bands is sensitive to the minimum 
energy rearrangement of the surface ions. Although this 
accounts for the observable responsiveness of TiO2(001), 
electronic charge-density analysis indicates a relatively 
lower reactivity with respect to defect surfaces. 

In interaction with nearest-neighbor oxygen, Pt derives 
s-p hybrid impurity states from the valence-bands, which 
lie above the VBM close to the Fermi level, as well as 
the ones from the conduction-bands. The character of 
these defect states controls the narrowing of the energy- 
band gap. In addition to the Pt-0 interaction, impurity- 
impurity coupling, depending on the Pt concentration, 
influences the position and the dispersion of defect states. 
Therefore, 1 ML cases exhibit strongly dispersed bands in 
the gap with a narrowing that refers to infrared activity 
and ultimately to metallization for Pt-Til substitution 
and for interstitial Pt implantation. Having the highest 
symmetry adsorption mode, Pt at the midpoint between 
two 02c's at 0.25 ML also exhibits a band gap narrowing 
of 1.42 eV corresponding to infrared region similar to low 
band-gap Pt/TiC>2(001)-(l x 1) cases. Pt concentration, 
therefore, is not the only factor that controls the band 
gap narrowing. Our calculations show that the number 
of impurity driven states that fall into the gap is propor- 
tional to the coordination of Pt with 01 and 02. This 
explains the occurrence of 1.42 eV gap narrowing for 0.25 
ML Pt concentration while all the remaining cases for 
(2x2) represent visible activity. Moreover, excess non- 
bonding charge due to the interaction of Pt with only 
one 01 in the adsorptional models in Figs.E^a) and[5jb) 
brings rather flat going bands. 

At 1 ML concentration under oxygen rich conditions Pt 
substituted Ti02 tends to be more stable. As O chemical 
potential gets slightly lower, the phase adopts a (2x2) 
surface periodicity. Further decreasing po toward Ti- 
rich conditions, the most stable phase, dominant over 
the allowed range of experimental conditions, happens 
to be adsorption on the surface for 0.25 ML coverage (b- 
2x2), which is semiconducting in the visible region with 
a band gap narrowing of 1.03 eV. Consequently, anatase 
TiO2(001) can be functionalized for visible activity by Pt 
incorporation. We have shown how it relates to impurity 
concentration, to coordination number of Pt with the 
near surface oxygen and to thermodynamic conditions 
which determine the chemical potential of constituent 
species, for useful applications. 
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